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Abstract 

We present a chaplygin gas Friedmann-Robertson- Walker quantum cosmological model. In this work 
the Schutz's variational formalism is applied with positive, negative, and zero constant spatial curvature. 
In this approach the notion of time can be recovered. These give rise to Schrodinger-Wheeler-DeWitt 
equation for the scale factor. We use the eigenfunctions in order to construct wave packets for each case. 
We study the time dependent behavior of the expectation value of the scale factor, using the many-worlds 
»' ! ' interpretations of quantum mechanics. 
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^ ■ 1 Introduction 
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In recent years supernova la (SNIa) observations show that the expansion of the universe is accelerating 
' [1] contrary to Friedmann-Robertson- Walker (FRW) cosmological models, with non-relativistic matter and 

radiation. Also cosmic microwave background radiation (CMBR) data [2, 3] is suggesting that the expansion 
of our universe seems to be in an accelerated state. This is referred to "dark energy" effect [41. Cosmological 
constant, A, as usual vacuum energy can be responsible for this evolution by providing a negative pressure 
[5, 6]. Unfortunately, the observed value of A is 120 orders of magnitude smaller than the one computed from 
field theory methods [5, 6]. Quintessence is an alternative to consider a dynamical vacuum energy [7], involving 
one or two scalar fields, some with potentials justified from supergravity theories [8]. However, the fine-tuning 
problem of these models which arise from cosmic coincidence issue has no satisfactory solution. 

The Chaplygin gas model is an interesting proposal [9], describing a transition from a universe filled with 
dust-like matter to an accelerating expanding stage. This model was later generalized in Ref. [10, 11]. The 
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generalized Chaplygin gas model is described by a perfect fluid obeying an exotic equation of state [11] 




where A is a positive constant and < a < 1. The standard Chaplygin gas [9] corresponds to a = 1. Some 
publications [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30] and reviews [24, 25] which 
studied the Chaplygin gas cosmological models have already appeared in the literature. 

Recently, quantum mechanical description of a FRW model with a generalized Chaplygin gas has been 
discussed in Ref. [31] in order to retrieve explicit mathematical expressions for the different quantum mechanical 
states and determine the transition probabilities towards an accelerated stage. In this paper we investigate the 
existence of singularities at quantum level in Chaplygin gas cosmological models. In the quantum cosmology 
the Wheeler-DeWitt (WD) equation in minisuperspace which determines the wave function of the Universe, 
can be constructed using ADM decomposition of the geometry [32] in the Hamiltonian formalism of general 
relativity. 

The presence of matter in quantum cosmology needs further consideration and can be described by funda- 
mental fields, as done in Ref. [33]. Using WKB approximation one can predict the behavior of the quantum 
universe which leads to determination of the trajectories in phase space. However, even in the minisuperspace, 
general exact solutions are hard to find, the Hilbert space structure is obscure and it is a subtle matter to 
recover the notion of a semiclassical time [34, 33]. 

In the present work, we describe matter as a chaplygin gas. This description is essentially semiclassical from 
the start, but it has the advantage of furnishing a variable, connected with the matter degrees of freedom, which 
can naturally be identified with time, leading to a well-defined Hilbert space structure. It is very convenient 
to construct a quantum chaplygin gas model. Schutz's formalism [35, 36] gives dynamics to the fluid degrees 
of freedom in interaction with the gravitational field. Using proper canonical transformations, at least one 
conjugate momentum operator associated with matter appears linearly in the action integral. Therefore, a 
Schrodinger-like equation can be obtained with the matter variable playing the role of time. 

Here, we use the formalism of quantum cosmology in order to quantize three Friedmann- Robertson- Walker 
chaplygin gas models in the presence of a negative cosmological constant. In Sec. 2 the quantum cosmological 
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model with a chaplygin gas as the matter content is constructed in Schutz's formalism [37], and the Schrodinger- 
Wheeler-DeWitt (SWD) equation in minisuperspace is written down to quantize the model. The wave- function 
depends on the scale factor a and on the canonical variable associated to the fluid, which in the Schutz 
variational formalism plays the role of time T. We separate the wave-function into two parts, one depending 
solely on the scale factor and the other depending only on the time. The solution in the time sector of the 
SWD equation is trivial, leading to imaginary exponentials of the type e~ lEt , where E is the system energy 
and t = —T. In Sec. 4 we construct wave packets from the eigenfunctions and compute the time-dependent 
expectation values of the scale factors. In Sec. 5, we present our conclusions. 



We need the Hamiltonian for a chaplygin gas model in the formalism developed by Schutz. The starting point 
is the action for gravity plus chaplygin gas, which in this formalism is written as 



where K ab is the extrinsic curvature and h a b is the induced metric over the three-dimensional spatial hyper- 
surface, which is the boundary dM of the four dimensional manifold M. Units are chosen such that the factor 
16ttG becomes equal to one. The first two terms were first obtained in [32]; the last term of (2) represents 
the matter contribution to the total action and p is the pressure. In Schutz's formalism [35, 36] the fluid's 
four-velocity is expressed in terms of five potentials <!>, £, f3, 9 and S 



where \x is the specific enthalpy. The variable S is the specific entropy, while the potentials ( and j3 are 
connected with rotation and are absent in models of the Fricdmann-Robertson- Walker (FRW) type. The 
variables $ and have no clear physical meaning. The four-velocity is subject to the normalization condition 



2 Model 




(2) 



u v = 



-{*, v + <;p, v + os iV ) 



(3) 



u v u v = — 1. 



(4) 



The FRW metric 



ds 2 = -N 2 (t)dt 2 + a 2 (t)g tj dx i dx j , 



(5) 
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is now inserted in the action (2). In this expression, N(t) is the lapse function and gij is the metric on the 
constant-curvature spatial section. 

Following the thermodynamic description of Ref . [38] , the basic thermodynamic relations take the form 



Po[l + H], h=l + U+p/p 



rdS = dH+pd(l/p ) 

(i + n)-< 



1 + a 



(i + n) 



1+a 



A 



Pl +a 



It then follows that to within a factor 



T = 



(i + n)-° 

1 + a 



s = (i + n)^ + -4^. 

Pa 



(6) 
(7) 

(8) 
(9) 



Therefore, the equation of state takes the form 

p = -A 
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S l/a 



The particle number density and energy density are, respectively, 



P = 



Pa = 



I i-^Z 

P + P 
h 



(10) 

(11) 
(12) 



where h = (6 + 6S)/N. After dropping the surface terms, the final reduced action takes the form 



S= I dt{ - 6— + 6kNa - Na 3 A 
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(13) 



The reduced action may be further simplified using canonical methods [38], resulting in the super-Hamiltonian 

2 l 

H = ~ ~ 6ka + + Aa3(1+a) ) 1+ " (I 4 ) 

where p a = —12da/N and — f§ • However, an analytical quantum mechanical treatment of this FRW 
minisupcrspace with the above Hamiltonian does not seem feasible. Therefore, it requires the following ap- 
proximation [31], 



1 Aa^ a+1 '> 1 1 
1 + ^- + 



1 + a Sp'+ a 2 1 + a V 1 + a J S2p2( 



- 1 



A 2 



2(l+a) 



6(a+l) 



(15) 



Hence, up to the leading order, the super-Hamiltonian takes the form 



V 2 

H= —£r -6ka + S— p* 
24a 



(16) 



The following additional canonical transformations, 



-(1 + a)p ip 1 S 1 +<xps, p T = S 1 +"p 9 , 



simplifies the super-Hamiltonian to, 



H=-^- a -6ka + PT ,. 



(17) 



(18) 



where the momentum pr is the only remaining canonical variable associated with matter. It appears linearly in 
the super-Hamiltonian. The parameter k defines the curvature of the spatial section, taking the values 0, 1, —1 
for a flat, positive-curvature or negative-curvature Universe, respectively. 

The classical dynamics is governed by the Hamilton equations, derived from Eq. (18) and Poisson brackets, 
namely 



' o= {a,NH} = -^, 



p a = { Pa ,NH} = -^p 2 a + 6Nk 



Via 

N 
24a 2 J 



(19) 



T = {T, NH} = N , 
p T = {pt,NH} = 0. 

We also have the constraint equation H = 0. Choosing the gauge N = 1, we have the following solutions for 
the system 

■- 2 k 



a 

2a ~ 2a' 







-6aa 2 — 6fca + pr- 



(20) 
(21) 



Imposing the standard quantization conditions on the canonical momenta and demanding that the super- 
Hamiltonian operator annihilate the wave function, we are led to the following SWD equation in minisuperspace 
(h = l) 

(22) 



— - - lUka 2 ^ + i2Aa— = . 
da 1 at 



In this equation, t = —T corresponds to the time coordinate. As discussed in [39, 40], in order for the 
Hamiltonian operator H to be self-adjoint the inner product of any two wave functions $ and must take the 



form 

poo 

($,*) = / a$*$>da, (23) 
Jo 

Moreover, the wave functions should satisfy the restrictive boundary conditions 

d#(a,i) 



*(0,t) = or 



= 0. (24) 

a=0 



da 

The SWD equation (22) can be solved by separation of variables as follows 

ip{a,t) = e- iEt ip(a) (25) 
where the a dependent part of the wave function (ip(a)) satisfies 

-ip"(a) + lUka 2 tp(a) = 24Eaip(a), (26) 
and the prime means derivative with respect to a. 

3 Results 

For k = the time-independent Wheeler-DeWitt equation (26) reduces to 

ip" + 24EatP = 0. (27) 
The Bcsscl functions are solutions of the above equation. Therefore, the time dependent solutions are as follows 



* £ = e - %Et VE 



A {V96E 3 

c\Ji\ — - — a 2 +C2Y1 — - — a 2 



(28) 



V 3 / 3\ 3 

Now, the wave packets can be constructed, by superposing these eigenf unctions with the following structure 

V{a,t)= A(E)V E (a,t)dE. (29) 
^0 

We choose c 2 = 0, for satisfying the first boundary condition (24). By choosing A{E) as a quasi-gaussian 
weight factor and defining r — V ' 9 ^ E , analytical expressions for the wavepacket can be found 

/■OO 

*(a,t) = y/E r v+1 e-~< r2+l £ r2t J„(rai)dr, (30) 
Jo 

where v = i and 7 is an arbitrary positive constant. The above integral is known [41], and the wave packet 
takes the form 



Figure 1 : The behavior of the expected value for the scale factor (a) (t) (solid line) and the classical scale factor 
a(t) (dashed line). 



where B = 7 — i-^t. Following the many worlds interpretation of quantum mechanics [42], we may write the 
expected value for the scale factor a as 

a^(a 1 t)*a^(a 1 t)da 



< a > (t) = 



J °° a^(a,t)*^(a,t)da 



which yields 



< a > (t) 



:t Z + 1 



(32) 



(33) 



_(32)V 

These solutions represent a bouncing Universe, with no singularity, which goes asymptotically to the corre- 
sponding flat classical models for late times (Fig. 1) 



(34) 



In the case k = 1 the time-independent Whcclcr-DeWitt equation (26) reduces to 



-ip"(a) + (-2AEa + 144a 2 ) V(a) = 0. 



(35) 



Defining new variable x = 12a — E we find 

dx 2 



x^_ 

144 + 144 



ip(a) = 0. 



(36) 



Equation (36) is formally identical to the time-independent Schrodingcr equation for a harmonic oscillator with 
unit mass and energy A 



.0+[-2A + W V]#r)=O, 



(37) 



where 2A = £? 2 /144 and w = 1/12. As much as the allowed values of A are n + 1/2, the possible values of E 
are 

E n = v/l2(2n+l) , rc = 0,l,2,... . (38) 

Thus, the stationary solutions are 

*„(a, t) = e~ iE " Vn (12« - £„) , (39) 

where 

Vn (x)=ir n (-^=)e- 2 /24 ) (40) 

with the n-th Hermite polynomial. The wave functions (39) are similar to the stationary quantum worm- 
holes as defined in [43]. However, neither of the boundary conditions (24) can be satisfied by the these wave 
functions. 

In the k = — 1 case the equation (26) reduces to 

ip"(a) + (2AEa + U4a 2 ) ip(a) = 0. (41) 

where the solutions are 

na,t) = e-^(12a + E)-^{c 1 M^ l{ ^±^) + C 2 W^ l( ^±^)j (42) 

48 '4 12 48 '4 YL J 

where M K ^\ and W K .\ are Whittaker functions. The Whittaker functions do not automatically vanish at a = 0. 
Therefore, in order to satisfy "5(0, t) = it is necessary to take both C\ ^ and Ci ^ 0, the same is applied 
to the second of the boundary conditions (24). 

4 Conclusions 

In this work we have investigated closed, flat, and open minisuperspace FRW quantum cosmological models 
(k = 1,0, —1) with chaplygin gas Universes. The use of Schutz's formalism for chaplygin gas allowed us to 
obtain a SWD equation in which the only remaining matter degree of freedom played the role of time. We 
have obtained eigenfunctions and therefore acceptable wave packets were constructed by appropriate linear 
combination of these eigenfunctions. The time evolution of the expectation value of the scale factor has been 
determined in the spirit of the many-worlds interpretation of quantum cosmology. 
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